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ABSTRACT 

Intensity scintillation in a laser beam at 0.63 

micrometers in the marine boundary layer has been studied 

over a 4.3 kilometer horizontal path across Monterey Bay 

and also from shore to ship at San Nicolas Island. Optical 

and microme teorological measurements of the refractive 

index structure constant, C^, agree to approximately one 

“ 8 

standard deviation giving values in the range 2.0 x 10 
to 5.0 x 10 (meters) ' . Optical and meteorological 

data. correlation improved as mean wind speed increased. 
Evidence was found to support the reported existence of a 
saturation region for the laser beam logarithmic amplitude 
.variance . 



2 



TABLE OF CONTENTS 



I. INTRODUCTION ------- 12 

II. PROJECT BACKGROUND ---------------- 15 

A. HISTORICAL SUMMARY ------- 15 

B. RECENT DEVELOPMENTS 17 

III. THEORY - - - - 20 

A. GENERAL REMARKS- ---------------20 

B. TURBULENCE THEORY- -------------- 22 

C. DERIVATION OF THE SCALAR WAVE EQUATION - - - - 25 

D. RYTOV'S METHOD ------- - - - 28 

E. APPLICATION OF TURBULENCE THEORY TO 

THE PROPAGATION EQUATION ----------- 31 

F. THEORETICAL ASSUMPTIONS 33 

G. PROBLEM APPLICATIONS ------------- 34 

H. METEOROLOGICAL TURBULENCE MEASUREMENTS - - - - 35 

IV. EXPERIMENTAL APPARATUS AND PROCEDURES 38 

A. THE SCINTILLOMETER SYSTEM 38 

B. THE THERMOSONDE SYSTEM ----- 43 

" C. SCINTILLOMETER PERFORMANCE VERIFICATION 

EXPERIMENTS- ----------------- 44 

1. Frequency Response Verification 
Experiments- ---------------45 

2. Sciptillation Spectral Analysis- ----- 46 

3. Laser Beam Intensity Profile Analysis- - - 47 

4. Laser Power Stability Analysis ------ 48 

V. DATA ANALYSIS- ------- ------ 51 

A. FIRST MONTEREY BAY EXPERIMENT- - 5 2 



3 



B. SAN NICOLAS ISLAND EXPERIMENT - 

C. SECOND MONTEREY BAY EXPERIMENT 

1. Data Taken on May 2, 1973- 

2. Data Taken on May 3, 1973- 

3. Data Correlation - - - - - 

VI. GENERAL RESULTS- - 

VII. FUTURE RESEARCH SUGGESTIONS- - - - 

VIII. CONCLUSIONS 

APPENDIX A - GRAPHS 

1 

APPENDIX B - FIGURES ---------- 

APPENDIX C - DATA TABLES -------- 

BIBLIOGRAPHY - - 

INITIAL DISTRIBUTION LIST 

FORM DD 1473 



56 

58 

58 

58 

60 

62 

64 

66 

67 



- 91 
■104 
-10 7 
■108 



4 



LIST OF GRAPHS 



I. First Monterey Bay Experiment Data- -------67 

A. Optical Data- ----------------67 

B. Meteorological Data -------------67 

II. First Monterey Bay Experiment Data - Scatter 

Plot- ----------------------68 

III. Refractive Index Structure Constant versus 
Time of Day San Nicolas Island Experiment 
Data- ----------------------69 

r ^ 

IV. Second Monterey Bay Experiment Data -------70 

A. Optical Data 2 May -------------70 

B. Meteorological Data 2 May- ---------70 

C. Optical Data 3 May 0630 - 1200- ------ 71 

D. Meteorological Data 3 May 0630 - 1200 - - - 71 

E. Optical Data 3 May 1200 - 1800- ------ 72 

F. Meteorological Data 3 May 1200 - 1800 - - - 72 

G. Optical Data 3 May 1800 - 2100- ------ 73 

H. Meteorological Data 3 May 1800 - 2100 - - - 73 

V. Scatter Plots of Data Taken during the 

Experiment of '2 May - -- -- -- -- -- -- --74 

A. Scatter Plot All Data- -----------74 

B. Scatter Plot of Data Taken when the Mean 

Wind Speed Was Less than 7 mph- -------75 

C. Scatter Plot of Data Taken when the Mean 

Wind Speed Was Greater than 7 mph ------ 76 

VI. Scatter Plots of Data Taken during the 

Experiment of 3 May - -- -*-----------77 

A. Scatter Plot of All Data- ---------- 77 

B. Scatter Plot of Data Taken when the Mean 

Wind Speed Was Less than 3 mph- ------- 7*8 



5 



- 79 



VII. 

VIII. 



C. Scatter Plot of Data Taken when the Mean 
Wind Speed Was between 3 mph and 6 mph- - - 

D. Scatter Plot of Data Taken when the Mean 
Wind Speed Was Greater than 6 mph - - - - - 

Scatter Plot of All Data from All Experiments - 

Logarithmic Intensity Scintillation Spectrum- - 



/ ** // 
3 7 



- 80 
- 81 
- 82 



6 



LIST OF FIGURES 



1. The Laser Scintillometer System- 

2. Laser Output and Pointing System ---------- 

3. Optical Detector Package -------------- 

4. The Laser Scintillometer -------------- 

5. Oscilloscope Trace of the Scintillation Signal - - - 

6. Multichannel Pulse Height Analyzer Display of the 
Probability Distribution of the Log- Intensity of 
the Laser Beam ------------------- 

,/ 

7. Scintilla cion Spectrum Analysis- - -- -- -- -- - 

8. Scintillation Signal Frequency Component Intensity 
versus Time- -------------------- 

9. Monterey Bay -------------------- 

10. San Nicolas Island ----------------- 



83 

84 

85 

86 

87 

87 

88 

89 

90 

90 



7 



LIST OF DATA TABLES 



I. First Monterey Bay Experiment Data --------91 

II. San Nicolas Island Experiment Data --------93 

III. Second Monterey Bay Experiment Data- -------94 



' r-i 

y 



/! 

fj 

J’ 



8 



ACKNOWLEDGEMENTS 



The scope of this project required the participation of 
a large number of personnel. The author would like to 
thank those people without whose support, this research would 
not have been possible. 

The project benefited from the support of the Office of 
Naval Research through the Naval Postgraduate School's 
Foundation Research Program and the Naval Ordnance Labora- 

S'* i 

tory. The cr operation of the managements of the Monterey 
Holiday Inn and the Hopkins Marine Station (of Stanford 
University)' was greatly appreciated. They provided rooms 
and electrical power for the laser and scintillometer 
systems . 

The Directorate of Industrial Operations at Fort Ord 
provided the radios used for establishing communications 
between measurement stations. Without these radios, coordi- 
nation of the measurements would have been nearly impossible. 

During the San Nicolas Island experiment, the assistance 
rendered by Mr. G. Gibbs, Mr. J. Karney and Mr. M. Oberg 
of the Naval Missile Center, Point Mugu, was invaluable. 

They arranged for range clearances, made independent optical 
measurements of atmospheric turbulence and assisted the 
optical research team in making ship- to -shore scintillation 
measurements . 

The Oceanographer of the Navy provided the facilities 
of the oceanographic research vessel Acania. The cooperation 



9 



of the Naval Postgraduate School's Oceanography Department 
and the captain and crew of the Acania was greatly 
appreciated . 

The meteorological research group was headed by 
Dr. Kenneth Davidson of the Meteorology Department. 

Included in this group were Dr. Thomas Houlihan of the 
Mechanical Engineering Department and Dr. Noel Boston of 
the Oceanography Department. This group has been responsi- 
ble for the meteorological experimentation and analysis 
constituting , he other half of the coordinated propagation 
study project, of which the work here described has been a 
part. The significance of the results is equally dependent 
on the meteorological and optical measurements. 

The optical research team was headed by Professor 
Eugene Crittendon of the Physics Department until his de- 
parture for Australia in January. Professor Alfred Cooper 
headed the optical team in his absence. Included in the 
optical research team were Professor Fred Buskirk, Professor 
Sydney Kalmbach and Professor Raymond Kelly, all of the 
Physics Department. These men spent many hours designing and 
running the optical experiments. In addition. Lieutenant 
Arthur Schroeder and Lieutenant David Beall, both 
graduate students in physics, gave many hours of their 
time helping to make measurements during the second 
Monterey Bay experiment. These men will be continuing the 
optical research effort as their masters thesis project. 

The author wishes them much success. 



10 



Technical help for the project was provided by 
Mr. Robert Moehler, Mr. Lynn May and Mr. Thomas Christian. 
Their practical- expertise in making systems work properly 
was greatly appreciated. 

A special word of thanks is richly deserved by my 
thesis advisor. Dr. Alfred Cooper. His participation, 
help and encouragement in the face of numerous "crises 1 ’ 
contributed a large part to the results and successes 
of this research effort. 

f a i 

Finally, ny wife, Valerie, deserves special acknowledge 
ment for her efforts in transporting equipment and 
especially for her moral support. 



11 



I. INTRODUCTION 



Atmospheric effects on laser beam propagation are 
divided into three general categories: refraction, absorp- 

tion, and particulate scattering. This investigation was 
exclusively concerned with atmospheric refraction effects. 
Refraction may be further subdivided into large-scale 
refraction, where the beam's spatial path is altered, and 
small-scale ref‘ v ac f :ion, where the beam is differentially 
focused and aefocused. Large-scale refraction, commonly 
referred to as beam wander, was found to be negligible. All 
research efforts were directed toward the small-scale refrac- 
tion phenomenon commonly referred to as beam or atmospheric 
scintillation. 

The major purpose of the research was to make and to 
correlate direct optical measurements of laser beam scintil- 
lation with meteorological measurements of certain atmospheric 
parameters. Using Rytov's method as described by Tatarski 
[1] , the optical and meteorological measurements of beam 
scintillation were compared by calculating the atmospheric 
refractive index structure constant which statistically 
describes atmospheric refractive index fluctuations. The 
research team attempted to predict laser beam scintillation 
by making micrometeorological measurements of atmospheric 
parameters . 

The project effort consisted of two primary areas, an 
optical determination of the refractive index structure 
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constant (C^q) and a meteorological determination of the 
refractive index structure constant (C^y) • These measure- 
ments were made simultaneously and compared for correlation. 
This experiment differs from previous experiments conducted 
by Livingston [2], Ochs [3] and Fitzmaurice [4], in that the 
laser path was over the ocean instead of over land. 

This research project was funded by the Naval Ordnance 
Laboratory and was divided into two related areas, meteorol- 
ogy and optics. _ The meteorological research group, headed 
' x 

by Dr. Kenneth Da/idson of the Naval Postgraduate School, 
made and evaluated the microme teorological measurements of 
atmospheric turbulence. The optical research group, 
headed by Dr. Eugene Crittendon (until his departure for 
Australia in January) and later by Dr. Alfred Cooper, made 
and evaluated the optical measurements of laser beam 
scintillation. 

Due to the wide extent of each of the areas of investiga- 
tion, each group worked independently of the other to obtain 
an estimate of. atmospheric turbulence. These estimates were 
compared for correlation after data reduction and analysis. 
The inherent nature of this organization promoted speciali- 
zation of the personnel involved in either the optical or 
meteorological areas of research. The author, by virtue of 
his background in electro-optics, concentrated his efforts 
in optical research. Hence, this report is devoted pri- 
marily to the procedures and conclusions of the optical 
research group. This thesis represents a comprehensive 
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description of the efforts and results of the optical 
propagation study. Meteorological techniques and experi- 
mental results are described briefly when they relate to 
optical data. 



II. PROJECT BACKGROUND 



A. HISTORICAL SUMMARY 

The earliest observations of scintillation in a light 
beam propagating through the atmosphere were made by ancient 
astronomers who noticed the "twinkling" of stellar images. 
Stellar scintillation was measured quantitatively by Middleton 
in 1949 [5] and 1952 [6]. Chandrasekhar, in 1952 [7], 
estimated the phase and amplitude fluctuations of starlight 
caused by atr. aspheric turbulence by means of a geometrical 
optics approach. In 1953, Mintzer [8] applied the Born 
approximation to the scalar wave equation to describe the 
propagation of sound waves through a turbulent ocean, a 
problem similar to the propagation of light through a turbu- 
lent atmosphere. 

Muchmore and Wheelon [9], in 1955, used Chandrasekhar’s 
approach to study line of sight propagation of radio waves 
in the troposphere. They calculated the variance and 
spatial correlation functions for the amplitude and the 
phase of the propagating wave. Two articles by Wheelon 
published in 1959 [10, 11] review the state of the art of 
wave propagation at that time. These articles cover 
turbulence theory' and atmospheric propagation of radio 
waves which is a similar problem to propagation of optical 
beams . 

In 1960 and 1961, the Russian monographs by Chernov [12] 
and Tatarski [1] were translated into English. These works 
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described the method proposed by Rytov [14] (see Section 
III.D for the derivation) for analyzing the diffraction of 
light by ultrasonic waves and later applied by Obukov [15] 
to the propagation of light waves through a turbulent 
atmosphere. Strohbehn [16] presents a good discussion on 
the range of validity and approximations made by Rytov' s 
method. The calculations in this project were based on 
Tatarski's formulation of Rytov' s method. Rytov' s method 
takes into consideration multiple scattering and diffrac- 

f'r f 

tion effects .nd thereby represents a significant advance- 
ment in the theory of atmospheric propagation of planar or 
spherical electromagnetic waves. Strohbehn' s article, 
published in 1968, is an excellent review of atmospheric 
propagation theory. 

Lumley and Panofsky [17] wrote (in 1964) a monograph 
covering many facets of turbulence theory as applied spe- 
cifically to the atmosphere. This defined the -meteorological 
theory used by the Naval Postgraduate School's research group 
to obtain estimates of the atmospheric refractive index struc- 
ture constant. In 1972, Davidson, in a Naval Postgraduate 
School report [18], comprehensively described the meteoro- 
logical measurements made and the theory applied in this 
project. Articles by Wyngaard (£t al . ) in 1971 and Panofsky [19, 
20] in 1968 provide further references on meteorological 
turbulence theory as applied by the research group. 
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B. RECENT DEVELOPMENTS 



In his monograph, Wave Propagation in a Turbulent 
Medium [1] , Tatarski reported the results of some of the 
earliest experiments that compared micrometeorological 
estimates to optical estimates of atmospheric turbulence 
using Rytov's method. These experiments were performed in 
the late 1950s using an incoherent light source. The 
propagation path ranged from 250 to 2000 meters over a 
Russian steppe. Experimental results strongly supported 
his formulation of propagation theory. The fluctuations of 
beam intensity and phase caused by atmospheric turbulence 
were found to have a log-normal probability distribution. 
Tatarski' s description of the scintillation frequency spectrum 
and the dependence of the fluctuation variance on path length 
were found to match propagation measurements. 

In 1966, Schmeltzer [22] published a derivation of the 
covariance of the logarithmic intensity and phase fluctua- . 
tions of a propagating beam of light. His results were 
more general than Tatarski's as they are valid for colli- 
mated and focused beams as well as plane and spherical 
waves. However, his results were in the form of compli- 
cated double integral equations and so were not directly 
useable. Fried evaluated these integrals to obtain analytic 
expressions for the logarithmic intensity variance of a spheri- 
cal wave [23] and for a collimated or focused beam [24]. His 
expression for spherical wave scintillation was in agreement 
with Tatarski's results. 
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In 1969, Ochs [3] repeated Tatarski’s experiments in 
Colorado using a helium-neon laser as a light source. His 
experiment found several deviations from Tatarski’s theory 
including the existence of a saturation region for the opti- 
cal measurements. When atmospheric turbulence was very 
high, Ochs found that meteorological measurements would 
continue to measure turbulence values but the optical 
measurements would saturate at a lower turbulence value. 

Optical and meteorological measurements would not corre- 

/ ; « 

late until tl j turbulence had decreased below the satura- 
tion level of the optical measurements. 

Similar experiments were performed over land by 
Fitzmaurice in 1970 [4] and by Dowling and Livingston 
in 1971 [2] . Both groups measured the scintillation of 
0.6328 and 10.6 micrometer laser beams and compared optical 
and micrometeorological estimates of atmospheric turbulence. ' 
Both groups found evidence of a saturation region. 

The work described in this thesis was part of an inter- 
departmental program formed as a combination of the efforts 
of existing optical and meteorological research groups at 
the Naval Postgraduate School. This study applied Tatarski’s 
formulation of Rytov’s method to describe atmospheric turbu- 
lence. Optical aiid micrometeorological measurements of 
laser beam scintillation and atmospheric temperature fluctu- 
ations have been made independently and then analyzed. Both 
types of measurement yielded values for the index of refraction 



18 



structure constant, C^. Theory predicts that if the 
turbulence is isotropic and homogeneous, these values of 
should be equal. However, two major unresolved questions 
remained. First, does saturation of optical values of 
C N exist in the high turbulence regime? If so, what param- 
eters determine this saturation? Ochs [3] reported that the 
saturation level seemed to vary with mean wind speed and 
the time of day. 

Second, Friehe, of Scripps Institute, reported recently 
that humidity fluctuations, which were neglected by Tatarski, 
may contribute significantly to refractive index fluctuations 
in the atmosphere [25], His findings have not been general- 
ly accepted by propagation researchers. The frequency response 
of available humidiometers is inadequate to check these pre- 
dictions . 

Although considerable evidence exists to support Tatarski* s 
propagation theory over land, the problem of propagation in 
the marine environment has not been resolved. If humidity 
variations influence refractive index fluctuations, laser 
propagation over the ocean may not be described by Tatarski* s 
formulation of Rytov*s method. This thesis examines the 
validity of Tatarski *s propagation theory in the marine 
boundary layer. 
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III. THEORY 



A. GENERAL REMARKS 

Atmospheric turbulence is a three-dimensional, nonlinear, 
stochastic process that occurs in air whenever certain fluid- 
mechanical parameters are exceeded. Random, dissipative 
fluctuations in the homogeneity of the atmosphere will cause 
the velocity profiles of the air to vary in time and space. 
These fluctuations have significant effects on the propaga- 
tion of a laser beam by continuously changing atmospheric 
density. This causes the refractive index to fluctuate 
with time and space. 

Atmospheric turbulence (linear) effects on laser beam 
propagation may be classified into two categories. Large- 
scale inhomogeneities will uniformly refract the laser beam 
and cause the beam to wander in space. Small-scale inhomo- 
geneities will differentially refract the laser beam and 
cause interference phenomena or scintillation. Small-scale 
turbulence alters the optical path length over different 
segments of the beam and generates interference patterns 
between adjacent wave fronts. This alters the initial phase 
and amplitude distributions of the propagating beam and 
destroys spatial and temporal coherence. 

Due to the inherent causes of turbulence, it is impossible 
to predict the transient velocity fields that arise for a 
given set of boundary flow conditions. Statistical analysis 
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must be used to describe atmospheric turbulence. By 

f 

analyzing the variations in atmospheric homogeneity, one 
may describe the fluctuations of the index of refraction 
over a given path. 

Relating the fluctuations in a propagating laser beam 
to atmospheric turbulence is a three-part problem. The 
first part is to find a convenient formulation for charac- 
terizing atmospheric turbulence. This may be done by 
using statistical estimators of the covariance and the 
structure function to describe fluctuations in the index 
of refraction. The second part of the problem is to derive 
the equations that govern the propagation of electromagnetic 
waves through a random medium, starting with Maxwell's 
equations. These equations will relate the magnitude and 
phase of the radiation to the index of refraction of the 
optical path. The final part of the problem is to combine 
the statistical description of the index of refraction with 
the propagation equation. The result will be an analytic 
relationship between the statistical description of the 
refractive index and a statistical description of the 
phase and amplitude of the propagating beam. 

The entire derivation is qxtremely long and difficult. 
Consequently only* the physically most important points of 
the derivation are discussed here, and much of the turbu- 
lence theory and mathematics are omitted. A complete 

treatment of the propagation problem may be found in 
references 1, 17, 22, 23 and 24. 
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In order to minimize the effects of beam wander on the 
scintillation measurements, the laser beam was diverged. 

The propagating diverged beam was approximated by a spheri- 
cal wave. Thus the general propagation equations obtained 
were evaluated for the case of a spherical wave. 

B. TURBULENCE THEORY 

\ 

Wave propagation in a random continuous medium is 

governed by a stochastic wave equation, whose index of 

refraction, n(r,t), is a continuous random point function, 

and characterizes the transmission medium. For line of 

sight propagation on the earth, optical paths are relatively 

short enabling one to say that the transit time of a given 

phase front between the source and receiver is much less 

than the time constant associated with changing spatial 

_ 3 

distribution of the refractive index (typically 10 seconds). 
The time dependence of n(r,t) will be accounted for by 
letting each n(r,t^), for i = 1, 2, 3... be considered a 
single realization of the random process, n(r). Hence, one 
can characterize the transmission medium by a refractive 
index in the form 

(1) n(?) = <n(r)) + n^r) = 1 + n-^r) 

where n^(r) represents the fluctuations in n(r). 

Tatarski [1] assumed that the power spectrum of the 
turbulence, $ n (K), can be described by the empirical 
Kolomogorov model where 
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( 2 ) 



* n (K) 



0.033 



C N 



K 



11/3 



exp( 



■ k 2 /k 



m 



C^j is a constant that parameterizes the total amount of 
energy in the turbulence. Tatarski showed that can be 
related to a statistical estimator by 

(3) D n (d) = C N 2 d 2/3 



where D (d) is the structure function of the refractive 
n 

index given by 



D n (d) = <[n (r) - n(r+d)] 2 > 



and hence 
•(4) 



2 - ([n(r) - n(r+d)] 2 / 
N ^773 



Assuming isotropic, homogeneous turbulence, the functional 
arguments r and d can be changed from vectors to scalar 
distances. 

The quantity K, the spatial wave number, is given by 
K = 2 tt /£ where l is commonly interpreted, in some sense, 
as the size of the turbulent eddy. This turbulence spectrum 
is valid for values of K that are in the inertial subrange, 
that is, the range where < % < Lq. The inner scale, 
characterizes the size of turbulent eddies where viscous 
dissipation of turbulent energy becomes predominant (for 
values of Z less than . K is the spatial wave number 
corresponding to the inner scale (i.e., K m = 2tt/£q) . Values 
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of £ greater than the outer scale, , characterize the 
region of the turbulence spectrum where the eddies are not 
homogeneous and where energy is added to the turbulence 
spectrum. For values of SL in the inertial subrange, the 
motion of the larger turbulent eddies causes the formation 
and transfer of energy to smaller eddies with very little 
viscous dissipation. These eddies then form and trans- 
fer energy to smaller eddies and so on until eventually, 
at approximately £ = £q, viscous effects become the pre- 
dominant method of energy dissipation and turbulent energy 
is converted to thermal energy. Typically, for atmospheric 
turbulence in the marine boundary layer, £^ is of the order 
of a few millimeters and Lq is of the order of a few meters 
(depending on the stability conditions and height above 
the boundary) . 

Turbulence theory shows that the covariance of the 
refractive index for homogeneous, stationary conditions 

(5) ’ C n (d) = <n(r+d) n(r)> 

is the three-dimensional Fourier transform of the turbulence 
energy spectrum 

(6) C (d)* = 4 it/ d /°°$ (K) K sin(Kd) dK . 

n o n 

Also, the structure function of the index of refraction, 
given by equations (3) and (4) , is related to the turbulence 
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energy spectrum by the transform relation 



(7) D n (d) = 8 it /” (I - Si K - d - Kd ~ ) -^n (K) 1(2 dK 

under conditions of isotropy, homogeneity and s tationarity . 
Given equations (2) and (6), one can relate C n (r) to C^. 

This means that can be used to describe the statistics 
of the fluctuations of the refractive index in the atmosphere. 



C. DERIVATION OF THE SCALAR WAVE EQUATION 

The equations describing the propagation of an electro- 
magnetic wave in the atmosphere are derived from Maxwell's 
equations for a non-conducting medium without charge sources. 



( 8 ) 



(a) 


V 


• D = 


0 


(b) 


V 


• B = 


0 


(c) 


V 


X E = 


3B 

at 


(d) 


V 


X H = 


3D 
3 1 


additional relations 


(e) . 


D 


= eE 


with e = 


(f) 


B 


= yH 


with y = 



constant . 



Taking the curl of equation (8c) one gets 

(9) V X V X E = V(V • E) - V 2 E = yV X || . 

Assuming a time dependence of e 1(i)t for all field quantities 
and substituting equation (8d) into the right hand side of 
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equation (9) , one finds 



(10) V(V • E) - V 2 E = -y X 



3 r 3D. 
3 t ^3 1^ 



= -ye 



3 2 E 

3t 2 



= yew E . 

Starting with Maxwell's equation (8a) and using a well known 
identify for the divergence operator, 

(11) V*D=V*eE=eV*E+E»Ve=0 
or 



( 12 ) 



V • E = 



E • Ve 



Substituting equation (12) into equation (10) gives 



(13) 



V 2 E + yew 2 E = -V( £ ) . 



The index of refraction for a medium is given by 



Cl 4) 



C , -.1/2 , r \ 1/2 

n = y c(ye) = (const) (e) 



where c equals the velocity of propagation in free space, 

V equals the velocity of propagation in the medium and the 
wave number is given by 



( 15 ) 



k = 

c 
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Hence 



(16) 



2 ,22 

co ye = (— y) (c ye) = k n 



and 



(17) 



V 2 E + k 2 n 2 E = -V( 6 * ■ ? e -) • 



Realizing from equation (14) that n = (const) e, then 

- 2 
rj ry ^ 

= -y- = 2Vln(n) and equation (17) becomes 



(18) V 2 E + k 2 n 2 E = -2V[E • Vln(n)] . 

The term -2V[E • Vln(n)] couples the various components 
of the field together and is responsible for polarization 
fluctuations in the linearly polarized beam. Fitzmaurice [4] 
shows that if the deviations of the refractive index from 
its mean value is very small (typically n-^/ (n> ~ 10 ^) and 
if the wavelength of the propagating wave is much less than 
the inner scale (£q) of the turbulence, then one can assume 
that the polarization fluctuation term is negligible. This 
gives the scalar wave equation 

(19) V 2 E + k 2 n 2 E = 0 
where E is any rectangular component- of E(r). 
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D. RYTOV'S METHOD 

At this point, Rytov [14] transformed equation (19) to 

(20) V 2 ln(E) + | Vln (E) | 2 + k 2 n 2 = 0 . ^ 

Now let 

(21) E (r) = A(r) exp[iS(r)] i/ 

and from equation (1), n(r) = 1 + n^ (r) . 

Define the complex phase, \p , as 

(22) E(r) = exp [ tp ( r ) ] = exp [lnA(r) + iS(r)] </ 
and therefore, 

(23) ln(E) = In (A) + iS = i|»(r). 

Substituting this result into equation (20) gives 

(24) V 2 rp + | Vip | 2 + k 2 (1 + n x ) 2 = 0. 

Now set 

(25) ^ = + so that 

(26) E = exp(i P Q + 4^) = E o exp(i[> 1 ) . 

E q is the unscattered wave which satisfies the equation 

(27) V 2 E q + k 2 E q = 0 / 
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or the transformed equation 



(28) V 2 ln(E Q ) + |Vln(E Q )| 2 + k 2 = 0. 

Substituting equation (22) into equation (28) , 

(29) V 2 iP o + | V^ q | 2 + k 2 = 0 . 

Now, substituting equation (25) into equation (24) gives 

(30) V 2 }p Q + V 2 i^ 1 + |vi» + V^ x | 2 + k 2 (l + 2 n 1 + n-j^ 2 ) = 0 

and subtracting this from equation (29) gives after some 
manipulation 

(31) V 2 ^ + j Vip x | 2 + 2 Vi i|» • Vip x + k 2 (2n^ + n^ 2 ) = 0. 

This is the nonlinear Ricatti equation. Since n^ is much 

2 2 2 

less than one, k n^ is much less than 2k n-^ and may be 
neglected. Similarly, since |V^^| << |V^ q | ~ k, then 
2 

| Vip x | is much less, than 2 V^ q • and one can also 

2 

neglect the . | | term. Finally, one can write equation 

(31) as 

(32) V 2 ip i + 2Vip o • Vijj 1 + 2k 2 n x = 0 . 

The assumption that | Vip | << |Vip 0 | ~ k is the only 
significant assumption made in Rytov's method. It is 
equivalent to saying that the change in the variation 
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of the complex phase, over a distance of about one 

wavelength in much less than 2ir radians. Realizing that 
the wavelength is much less than the inner scale, it 

is difficult to imagine the phase changing significantly 
in a propagation distance less than &q. This approximation 
is generally referred to as the Rytov approximation in the 
literature . 



Substituting equation (34) into equation (32) and expanding 
terms gives 



A new quantity is now defined 



/ 



(33) 



ip 

W = iKe 0 = rp.E 
1 r l o 



so that 



(34) 




(35) 



V 2 W + k 2 W = - 2k 2 n. E 



1 o * 



This differential equation has the solution 



(36) W(r) = f 2k 2 n. (r ' ) E (r‘) G(r - r') dV‘ 
v' 1 0 



where G(r - r') is the Green’s function 



(37) G(r - f') = exp[ik(? - ?')] 



for a point source. 



4tt(? - r') 
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Therefore , 



W 

(38) ^(r) = | 



o 



= 1 r 2 k 2 n-^r*) E 0 (r*) G(r- r') dV 




Equation (38) with n^(r) as a random continuous function 
characterizes the scintillation in the propagating wave if 
n-^(r) is known analytically. ^ 

It is not practical to measure the index of refraction 
at every point over the optical propagation path. Instead 
it is necessary to statistically characterize n-^(r) using 
the covariance function, C n (r). This brings one to the 
third part of the problem. The statistics of if^(r) must 
be calculated starting from equation (38) and the statisti- 
cal description of the covariance of the refractive index 
as given by equation (6). 

E. APPLICATION OF TURBULENCE THEORY TO THE PROPAGATION 
EQUATION 

After making many transformations, substitutions and 
integrations, Schmeltzer [22] obtained the statistical 

? 

expression for the logarithmic amplitude covariance, C^(d) 
for a beam propagating in a random medium which Fried [23] 
reduced to 



( 39 ) 





Z oo 

Re { / ds / $ (K 1/2 ) J (K 1/2 d -) 
o O ft 0 Z 



X [1 - exp ( 
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for the case of a spherical wave. The quantity "z" is 
the optical propagation path length and "s" is a dummy 
variable of integration with the dimensions of length. 

Fried integrates this expression for a spherical wave 
over s and K to obtain the relation 

I 

(40) C^ A (d=0) = C £ 2 = 0.124 k 7/6 z 11/6 C N 2 where 

(41) C M 2 = ([In (A) - In (A)] 2 ) 

and A is the amplitude of the wave function, |E(r)|. 

The logarithmic intensity variance, C {T 2, is a more readily 
measurable quantity 

(42) C £I 2 = ([In (I) - ln<I> ) 2 > 

2 

where I is the intensity of the radiation (I = A ). Thus, 

> 

(43) C £I 2 = ([ 21n (A) - 21n <A> ] 2 ) 

= 4([ln(A) - ln<&] 2 ) 

= 4 C £A 2 and 

(44) C u 2 = 4 [0 . 124 k 7/6 z 11/6 C N 2 ] . 
or finally 

(45) C N = 1.42 k" 7/12 z" 11/2 Cjj-j . . 

The foregoing is a brief summary of the major steps in 
the development of the governing equation for the propagation 
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statistics of a spherical wave. Much of the detailed 
mathematics involved in the derivation has been omitted. 

For the complete derivation, the original papers by 
Tatarski [1], Schmeltzer [22], and Fried [23] should be 
consulted. However, the approximations involved in 
reducing the equations to analytic form are interesting 
and will be analyzed here. 

F. THEORETICAL ASSUMPTIONS 

The major assumption that underlies the entire theory 
is that the turbulence energy spectrum can be described by 
the empirical Kolomogorov model [equation (4]]. This model 
was developed in connection with wind tunnel experiments. 
Meteorological measurements by Tatarski [1], by Wyngaard, 

Izumi and Collins [19] and by Panofsky [20] verified that 
it was applicable in the atmospheric boundary layer over 
the inertial subrange. -1 It was also implicitly assumed that 
only those turbulent eddies in the inertial subrange contribute 
significantly to beam scintillation. Eddies with l > Lq would 
tend to cause beam wander. Eddies with i < Lq have little 
effect on the scintillation because viscous forces cause 
rapid collapse. 

The assumption that X << was used to justify ignoring 
the cross polarization term in the scalar wave equation 
(equation 18) . The helium-neon laser used in the experiment 
had a wavelength of 0.6328 micrometers which is much less 
than the dimensions of the inner scale (£q was of the order 
of a few millimeters). 
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In order to use the perturbation technique to express 
the variations in the index of refraction, the fluctuations 
in refractive index were assumed to be much smaller than the 
mean value. Meteorological measurements confirm that the 
index of refraction in the atmosphere does not increase 
by more than a few per cent due to density variations; 
hence, this assumption seems justified. 

The assumption that z/ 1 ^ <<Y is a very weak restriction 

at optical wavelengths. For a helium-neon laser and an 
inner scale value of one millimeter, the propagation path 

g 

would have to be greater than approximately 10 meters for 
this assumption to become invalid. 

The final assumption was that turbulence was isotropic 
and homogeneous. This approximation is valid for large 
areas with a smooth regular topography. Ocean conditions 
are conducive to isotropic homogeneous turbulence. This 
approximation means that a point turbulence measurement can 
characterize the turbulence over an entire propagation path. 

The experimental .conditions seemed to satisfy the require 
ments of Tatarski's, Schmelt zer ' s and Fried's theory for the 
propagation of a spherical wave in a random medium. A goal 
of this experiment was to test the validity of their propa- 
tion theory, in the marine boundary layer. 

G. PROBLEM APPLICATIONS 

The calculation of from equation (45) required values 
of the wavelength of the laser beam, the optical path length 
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and the variance of the logarithmic intensity fluctuations 
of the laser beam. The wavelength of the laser was known 
to be 0.6328 micrometers. The optical propagation path 
length was measured on a nautical chart. Log- intensity 
variance was measured using the laser scintillometer which 
is described in the equipment section of this thesis. The 
calculated values of were designated C^q, indicating that 
they were derived from optical measurements. These values of 
Cno were compared with values of that were calculated from 
meteorological measurements (designated C^ T ) . A brief 
explanation of the theory and equations that describe the 
calculation of is appropriate at this point. The 

meteorological theory is described in detail by Panofsky 
•[20] and by Wyngaard, Izumi and Collins [19]. 

H. METEOROLOGICAL TURBULENCE MEASUREMENTS 

Variations in the index of refraction in the atmosphere 
are caused by variations in the density of the air. These 
variations are due to fluctuations in the temperature, 
pressure, humidity and perhaps gravity waves. Empirical 
measurements [4] show that the index of refraction at opti- 
cal wavelengths can be represented by' the equation 

(46) n (r) = 1.0 + 77.6 (1 + - ° -i ° - ° - 7 - 5 ) X 10 -6 

X 

where atmospheric pressure., P(r) is in millibars, absolute 
temperature, T(r), is in Kelvins, and wavelength, X, is in 



35 



micrometers. Evaluation of equation (46) at a wavelength 
of 0.6328 micrometers gives 



(47) 



n(r) = 1.0 + 79.0 X 10" 6 



It is generally accepted in atmospheric propagation 

i 

studies that temperature fluctuations are the predominant 

mechanism causing index of refraction fluctuations in the 

1 2 
atmosphere. Equation (4) define S C N as 



(4) 






2 , (fn(r) - nfr+d)l 2 ) 



TUT 



One can also define a temperature structure constant, C^, by 

r 2 , ( f T (r) - T(r*d^l 2 ) 



(48) 



iF? 



Rearranging terms gives 

P 



(50) 2 _ 

L N 



([(79. 



0 



T(r) 



X 10" 6 ) - (79.0 P 



T (r+D) 



X 10 



_6 )] 2 ) 



,2/3 



and assuming that the temperature variations are much smaller 
than the mean temperature, then 



(51) 



1 1 ~ 1 ~ 1 

TUT " T(r+'d')' ' TTrJ " T 



and 



^Recently, Friehe of Scripps Institute, in a private 
communication, suggested that humidity variations may contrib- 
ute significantly to refractive index fluctuations under con- 
ditions of high humidity. This has not been generally accepted. 
It is difficult to measure humidity fluctuations due to a lack 
of humidiometers that will measure fluctuations with a fre- 
quency greater than 15 Hz [25]. 
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(52) C 



((79.0 X 10' 6 ) 2 [T (r) - T(r+d)] 2 \ 



N 



TUT 



Substituting equation (48) into equation (52) gives 



(53) 



V = ( 79 -° ^2 x 10 ’ 6 ) 2 (C t 2 ) . 



Thus the analytic expression for is 



(54) 



'NT 



= (79.0 ^ X 10" 6 ) C T 



In order to calculate C^ T using equation (54) measurements 
of mean pressure, mean temperature and C^, as given in equation 
(48) must be made. These measurements will be described in 
the experiment section of this thesis. 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 



This section will describe the apparatus and procedures 
used to make measurements during the project. The specific 
results and conclusions derived from the measurements will 
be reported in later sections. 

Two major systems were used to measure the data obtained 
during the three main experiments. The scintillometer 
system measured the standard deviation of the logarithmic 
intensity distribution, C^j, for the transmitted laser beam. 
These measurements furnished the data for the calculation of 
C^jq from equation (45). The thermosonde system measured the 
temperature structure constant, C^, mean temperature, T, 
and the mean pressure, P. These measurements furnished 
the data used in the calculation of from equation (54). 

Additional experiments were performed to measure the intensity 
stability of the laser, the power spectrum of the scintilla- 
tion and the degree of truncation of the laser beam by the 
laser output optics. These experiments verified the 
scintillometer system's performance and accuracy. The 
equipment used in the experiments will be subsequently 
described in this section. 

A. THE SCINTILLOMETER SYSTEM 

The scintillometer system measured the logarithmic 
intensity variance of the transmitted laser beam. The 
main components of the system are illustrated in Figure 1. 
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A low powered helium-neon laser transmitted a beam to an 
optical receiver. The receiver detected the incoming beam 
and generated an electrical voltage proportional to the 
instantaneous intensity of the beam. This signal was ana- 
lyzed by the laser scintillometer and recorded on magnetic 
tape. The standard deviation of the logarithmic intensity 
of the laser beam was computed by analogue processing methods. 
This section analyzes in detail the operation of this system. 

The laser output and pointing system controlled the spot 
size and direction of the beam (Figure 2). The laser beam 
was focused by the first lens to a position near the focal 
point of the output lens. An adjustable mount was pro- 
vided to give the output lens two degrees of translational 
freedom. The spot size was controlled by changing the focal 
point of the output lens with respect to the fixed focal 
point of the first lens. During normal operation, the beam 
was diverged to minimize beam wander effects. A telescope 
with cross hairs, bore-sighted with the laser, provided an 
accurate way of aiming the beam. Beam alignment with the 
target was controlled by elevation and azimuthal adjustment 
screws. A chopper was provided to interrupt the beam for 
ten seconds of every minute during which adjustments were 
made on the background compensator at the receiver. 

After transmission through the atmosphere, a portion of 
the scintillating beam plus background radiation was collected 
by the input optics in the detector package (Figure 3) . The 
incoming radiation was passed through a 0.63 micrometer 
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bandpass filter and focused by a lens through two small 
apertures. The detector package was aligned in such a 
manner as to pass the laser beam plus background radiation 
through the first aperture to the first silicon photocell. 

The second aperture passed only the background signal which 
was reflected from a small prism to the second photocell. 

The photocells generated an output voltage proportional to 
the intensity of the radiation that impinged upon the photo- 
sensitive surfaces. 

The two signals were then subtracted in a differential 
amplifier (Figure 4) thus cancelling the background signal 
and passing only a signal that was linearly proportional to 
the laser radiation. During operation, the laser beam was 
interrupted for ten seconds every minute, during which time, 
the background signals from the two photocells could be 
electronically adjusted in order to exactly cancel each other. 

The scintillation signal (Figure 5) was amplified by a 
factor of five and then passed through a logarithmic ampli- 
fier. The logarithmic amplifier computed the logarithm of 
the original signal and sent it to a wide band power ampli- 
fier. The signal level was multiplied by a factor of ten 
and recorded on magnetic tape with an F.M. conversion 
instrumentation tape recorder. The signal was then fed to 
a multi-channel pulse height analyzer, which sampled the signal 
level periodically. The pulse amplitude was used as a 
criterion to address one set of 200 analyzer memory channels. 
These channels counted the number of pulses of a given 
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voltage level that were accumulated over a variable time 
period (usually around four minutes). The pulse height 
analyzer output consisted of a CRT display of the number of 
pulses of a given signal level as a function of the logarithm 
of the beam intensity. This corresponded directly to a non- 
normalized, probability distribution of the logarithm of 
the laser beam intensity, ln(I). The CRT display was 
photographed using a polaroid camera and the photographs 
(example shown in Figure 6) were later analyzed to determine 

C N0‘ 

To aid in the analysis, a calibration signal generator 
provided direct current calibration signals at 25, 50, 

100 and 200 millivolt levels which were accumulated and 
.stored on a second set of 200 pulse height analyzer channels. 
The two sets of 200 channels could be displayed in such 
a manner as to superimpose the calibration signals and the 
' scintillation pulse height distribution. Amplifier and 
pulse height analyzer memory drift were found to be negli- 
gible for the system. 

Measurements by Tatarski [1], Ochs [3] and Dowling [4] 
have shown that the probability distribution of the logarithm 
of beam intensity is Gaussian. This was expected as a 
direct result of the central limit theorem and the randomness 
of atmospheric turbulence. Realizing that the pulse height 
distribution was Gaussian, then the maximum point corresponded 
to the mean value of the distribution and the standard devia- 
tion of the logarithmic intensity, C^j, could be measured 
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directly from the photographs. Realizing that for a 
Gaussian distribution, the full width at half maximum, Xpi/2, 
was related to the standard deviation by 



(56) 



C &1 = 0.425 X pl / 2 



then the standard deviation could be determined from the 
following equation. 



"S" is a scaling parameter that related the distance along 
the abcissa in centimeters to the signal value in log- 
intensity units. S was determined by measuring the distance 
in centimeters between calibration marks and by applying 
the relation 



where V is the calibration voltage, X is the distance along 
the abcissa in centimeters and AX is the distance in centi- 
meters between two calibration marks, Vp and V^. Once C^j 



was determined, C^j was computed using equation (45) . 

Implicit in the system computation of was the 
approximation that ln([l) = (ln(I)^> . This approximation 



(57) 



C$,1 = °* 425 S ( x pl/2 in cm ) 



(58) 



[ln(V 2 ) - In (Vp) .] 



(X 2 cm - XpCm) 




(AX) cm 
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was tested by computing both quantities for typical 
probability distributions of intensity and comparing the 
difference. The error introduced for large values of C^q 
was less than 2% and was correspondingly smaller for lower 
values of C^q. This error was judged negligible compared 
to the experimental error involved in measuring C^, and 
was ignored. 

The dynamic range and frequency response of the 
scintillometer were wide enough to process all scintilla- 
tion signals without distortion. The logarithmic amplifier 
had a dynamic range of 80 db and distorted voltage levels 
from the differential amplifier that were less than one 
millivolt or greater than ten volts. By choosing an 
appropriate gain setting, the scintillation signal's mean 
value at the differential amplifier output terminal was 
controlled to be approximately 100 mV. All other system 
components had dynamic ranges greater than 80 db. The 
minimum signal levels detected during the experiments were 
around five millivolts while the maximum levels were below 
one volt. System frequency response was flat from dc to 
1.5 KHz. However, the scintillation power spectrum was 
measured to be negligible beyond 300 Hz and hence, the 
signal frequency response was undistorted. 

B. THE THERMOSONDE SYSTEM 

The meteorological equipment used for measuring 
consisted of a platinum wire thermosonde unit, a barometer, 
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a thermometer and a psychrometer . The platinum wire 
thermosonde unit measured the temperature structure function, 
C T , by measuring the changes in the difference of resistance 
between two exposed platinum wires mounted on probes. The 
resistance of each wire was proportional to the temperature 
of the wire. The system electronics used analogue processing 
techniques to compute C T from the equation 

c 2 _ ( f T(r) - TQH-d)] 2 ) 

T 

where the value of "d 1 , the probe separation distance, was 

constant. This system had a linear frequency response from 

dc to 1 KHz and a wide dynamic range. The system generated 

2 

a voltage signal corresponding to which was fed to a 

strip chart recorder. The thermosonde output was recorded 
as a function of time. 

Mean pressure was recorded as a function of time from 
a barometer. Mean temperature was read and recorded from a 
mercury thermometer. Relative humidity was also read and 
recorded from a psychrometer. was calculated from 

measurements of C^,, mean temperature and mean pressure 
using equation (54) . 

C. SCINTILLOMETER PERFORMANCE VERIFICATION EXPERIMENTS 
Three minor experiments were performed in conjunction 
with the scintillation measurements to determine the 
scintillation power spectrum, laser output intensity 
stability and beam truncation by the output optics. 
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Each experiment was necessary to verify the performance of 
the scintillometer system. 

1 . Frequency Response and Dynamic Range 

In order to assure that all laser beam signal levels 
were being accumulated by the pulse height analyzer, the 
dynamic range and frequency response of the system were 
tested. The system response to a wide range of input signal 
levels was found to be logarithmic for signals greater than 
one millivolt and less than ten volts. Typical log-intensity 
distributions ranged from approximately 15 to 400 millivolts 
and all measured distributions were within the system's 
dynamic range. 

Similarly, frequency response tests showed that the 
system would record and process signals with a frequency 
range of dc to 1.5 KHz without distortion. It was then 
necessary to determine the scintillation power spectrum to 
justify the frequency response of the scintillometer system. 
An analogue Fourier spectrum analyzer was designed and 
assembled using a wave analyzer, calibration oscillator and 
an X-Y plotter (Figure 7) . 

The FM instrumentation tape recorder fed a measured 
scintillation signal through a power amplifier to a wave 
analyzer. The wave analyzer passed the signal through a 
tunable bandpass filter (bandwidth of 3 Hz) and then 
generated a voltage equal to the rms value of the filtered 
signal. The rms voltage signal was recorded as a function 



45 



of time on an X-Y plotter. The system output was a graph 
of the scintillation signal instantaneous relative log- 
intensity at the specific frequency as a function of time 
(Figure 8) . The graph was then mechanically integrated 
(using a planimeter) to find the average relative intensity 
of the signal at the specific frequency. By tuning the 
bandpass filter to a variety of frequencies and recording 
the signal intensity, the scintillation intensity power 
spectrum could be measured. 

The tunable bandpass filter had a low frequency 
cut-off at 20 Hz. The system frequency response was ex- 
tended to 5 Hz by playing the tape of the scintillation 
signal at four times the original recording speed and thereby 
increasing the original signal frequency by a factor of four. 

A tunable oscillator was used to calibrate the wave 
analyzer bandpass filter by generating a pure sinusoidal 
wave at the desired frequency and tuning the filter to match 
that frequency. An oscilloscope was used to measure the 
oscillator's output accurately. 

2. Scintillation Spectral Analysis 

The scintillation power spectrum was analyzed to 
verify that the system did not distort the high frequency 
components of the scintillation signal. Spectrum analysis 
of three one-minute duration signals was accomplished using 
the Fourier spectrum analysis system. The three signals 
analyzed represented the maximum, minimum and average values 
of recorded during the first Monterey Bay experiment-. 
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The relative intensities of the various frequency components 
were averaged to give a representative scintillation spectrum 
as shown in Graph VIII. 

High frequency components above approximately 300 Hz 
had negligible power levels. The system's high frequency 
distortion limit was approximately 1500 Hz. This was 
determined by the tape recorder frequency response. The 
experiment showed that the high frequency components of the 
scintillation signal were not distorted by the system 
electronics . 

3. Laser Beam Intensity Profile Analysis 

The intensity profile of the laser beam was examined 
to determine the degree of truncation by the output optics. 

A scanning photomultiplier was used to measure beam intensity 
as a function of position in the plane of the output aperture. 
Intensity profiles were measured with and without the output 
lens in position and examined for beam truncation by the 
lens holders. 

The intensity scans showed this profile to be 
essentially Gaussian with negligible truncation. The pro- 
files also indicated that the laser power level varied in 
time, during the warm-up period. Further analysis of 
intensity stability was necessary to show that the variations 
in the radiation received by the scintillometer were due to 
atmospheric propagation effects and not to internal intensity 
variations . 
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4 . Laser Power Stability Analysis 

Laser intensity variations were monitored and 
recorded using a photomultiplier tube and strip chart 
recorder. Measurements of intensity as a function of time 
were taken in the center and wings of the beam profile. 

All measurements showed large fluctuations in beam intensity 
during the warm-up period. These fluctuations were in the 
form of a periodic modulation that gradually decreased in 
amplitude and increased in period. Maximum modulation 
occurred immediately after lasing began. After approxi- 
mately one hour of operation, the intensity stabilized and 
the modulation became negligible. However, it was noticed 
that a 5% modulation could be induced in the beam by blowing 
air across the laser vents after the laser had reached 
thermal equilibrium. The cool air created thermal stresses 
on the laser tube and mirror supports, and caused variations 
in the optical gain. 

Realizing that during the second Monterey Bay experi- 
ment, the laser's operating environment continuously subjected 
the laser to thermal stress from meteorological phenomena, 
the validity of the data would be in doubt if the internal 
modulations were large. To test the effects of an outdoor 
environment on thb intensity stability, the laser, photo- 
multiplier and strip chart recorder Were assembled on the 
roof of Spanagel Hall and the experiment was repeated. A 
short propagation path minimized atmospheric scintillation 
in the beam. 
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The laser intensity fluctuations stabilized from 351 
modulation when lasing commenced, to approximately 4% 
modulation after one hour. The level of modulation varied 
with wind velocity but never exceeded 41. 

Assuming that a 5% modulation was the maximum 
internal intensity variation encountered during the second 
Monterey Bay experiment, the magnitude of the perturbation 
effect on the measured value of C^q was calculated from 

(59) C u 2 = {[ln(I/(l) )] 2 j where 

I = I + AI 
s m 

In the equation, I g is the instantaneous intensity of the 
received beam due to scintillation, and AI was the devia- 
tion from the average intensity due to internal perturba- 
tions of the laser output intensity. Calculations using 
equation (59) with values of AI^ = ±0.05^1^ and a typical 
probability distribution of I g showed that the maximum 

increase in C^q values was 6.41 (computed for a of 

“8 “1/3 

2.0 x 10 m 7 , calculations showed that the perturbation 

effect would be less than 2.221. 

. - • These calculations represented the maximum increase 
in C^o values due to a continuous, random perturbation of 
AI m = ±0.05 (i) . In actual conditions, the measured average 
intensity perturbation was less than 5%. These estimates 
are an upper bound for possible perturbation effect errors 
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in the data taken on the second and third of May. The 
laser was shielded inside a motel room during the first 
Monterey Bay experiment and perturbation effects should have 
been negligible. Perturbation effects on the laser used 
during the San Nicolas Island experiment by the research 
group from the Naval Missile Center are an unknown quantity. 
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V. DATA ANALYSIS 



The primary goals of the laser propagation research 
project were: 

(1) to develop instrumentation and expertise for 
measurement of laser beam scintillation and micro- 
meteorological turbulence, 

(2) to examine the correlation between micrometeoro- 
logical atmospheric turbulence and laser beam 
scintillation and 

(3) to verify Tatarski's [1] formulation of the theory 
of laser propagation through a turbulent atmosphere, 
in the marine boundary layer. 

Three main experiments were performed between January 
and May of 1973 during which optical and meteorological 
data were collected. The first experiment was conducted 
on January 10. The laser radiation was transmitted across 
Monterey Bay from the Holiday Inn to. a receiver in a labora- 
tory at Hopkins Marine Station (see Figure 9) . Micrometeoro- 
logical measurements were made on-board the Naval Postgraduate 
School's Research Vessel Acania . The second experiment was 
conducted at San Nicolas Island on March 3. A laser on the 
Island was aimed at an optical detector mounted on the Acania . 
An additional set of laser scintillation measurements was 
made by a research group from the Naval Missile Center, 

Point Mugu, across Laser Bay on San Nicolas Island (see 
Figure 10) . The third experiment was conducted across 
Monterey Bay between the Holiday Inn and Hopkins Marine 

t 

Station on the second and third of May. Micrometeorological 
measurements were taken on the beach adjacent to the 
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Holiday Inn. The data collected during the three main 
experiments will be reported and analyzed in this section. 
General conclusions based on the results of the three 
experiments will be presented in the following section. 

A. FIRST MONTEREY BAY EXPERIMENT 

The first Monterey Bay experiment was carried out on 

January 10. The laser output and pointing system was 

) 

placed inside a room in the Monterey Holiday Inn. The beam 
was transmitted a distance of 4.3 kilometers across the Bay 
to the receiver located at Hopkins Marine Station. The 
transmission path was horizontal and approximately six meters 
above the surface. This path was closely transverse to the 
prevailing wind direction, and was chosen for this reason 
and to minimize land effects on turbulence homogeneity. The 
meteorological data were taken on-board the R.V. Acania 
which was anchored at approximately the mid point of the 
laser transmission path. The optical and micrometeorological 
measurements were- taken at the same height above the sea 
surface . 

The optical and meteorological estimates of are 
listed in Data Table I and are shown as functions of time on 
Graphs I A and B. Graph II is a scatter plot of the values 
of C NQ plotted against C NT> 

The optical measurements of were made in the 
laboratory by processing the tape recorded scintillation 
signal. The pulse height analyzer accumulated the signal 
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probability distribution over consecutive five-minute 
periods. This corresponded to a five-minute integration of 
the laser beam scintillation signal. Meteorological data 
were recorded on magnetic tape and later transferred to a 
strip chart from which the average values of were 
measured for five-minute periods. Lack of experience in 
operating the scintillometer system caused several inter- 
ruptions in data collection. 

The experimental error in the optical measurements of 

Cnq was predominantly due to errors involved in measuring the 

full width of the pulse height distribution. For values of 

-8 “ 1/3 

Cj^q between 2 and 5 (X10 m ' ) , the standard deviation of 
the measurements was approximately 5%. Instrumentation 
error and the error involved in measuring optical wavelength 
and optical path length were much less than 5%. 

The experimental error involved in the meteorological 
measurements of was predominantly due to inaccuracies in 
determining the average value of over five-minute intervals 
from strip chart recordings of the instantaneous measurements 
of Cj. To complicate this process, the thermosonde unit 
had high and low scales for different extremes of turbulence. 
Whenever the scale changed, 30 seconds of data were lost 
due to system transients. The standard deviation was 
estimated as approximately 20% for the measurements. 

Errors involved in measuring pressure and temperature and 
internal system errors were negligible compared to averaging 
errors. Experimental error bars are indicated at one point 
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on each graph. These estimates of experimental error 
apply to all points on the graph. 

In analyzing the data, it is important to realize the 
fundamental difference between the nature of the meteorologi- 
cal and optical measurements. The meteorological measure- 
ments represent the degree of turbulence at only one point 
along the propagation path. The optical measurements 
represent the integration of turbulence effects on laser 
propagation over the total propagation path. The point 
measurements of turbulence showed greater variations in the 
values of because localized cells of high turbulence 

strongly affected the meteorological sensors. The optical 

\ 

measurements did not respond to localized cells of strong 
turbulence unless the dimensions of those cells were sig- 
nificant compared to the propagation path length. The laser 
beam scintillation, in effect, averaged the effects of 
atmospheric turbulence over the entire optical path. This 
is equivalent to saying that the assumption of homogeneous 
turbulence is not entirely applicable to distances on the 
order of several kilometers. Hence, one should analyze the 
data by comparing the average level of C^ T to the optical 

measurements of C xiri . 

NO 

The data showed a strong correlation between the optical 

and meteorological values of with' the exception of three 

data regions on Graph I. At approximately 1325, the value 

- 8 - 1 / ^ 

of increased to 12.0 X 10 meters but the corre- 

- 8 _ i / 

sponding value of C^q increased to only 4.36 x 10 meters ' 
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increased 



Again at approximately 1350, the value of 
-8 “1/3 

to 15.0 x 10 m ' whereas C^q remained steady at 
3.25 x 10 ^m Finally, between 1525 and 1630, the 

values of showed large variations. Corresponding optical 

data remained relatively constant. The average value of 
the meteorological data was approximately equal to the 
optical data. 

These deviations between and may have been 

caused by inhomogeneous turbulence which strongly affected 
the temperature sensors but was not sufficiently wide- 
spread to affect the laser beam propagation. The fact that 
the average value of (between 1525 and 1630) was approxi- 

mately equal to the relatively steady value of C^q tended 
to support this explanation. However, it was also possible 
that the turbulence was homogeneous and that the optical 
values of became saturated. Although the saturation 
effect has been observed in over- land measurements by Ochs 
in 1969 [3] , the cause of saturation has not been clearly 
understood or explained. 

Graph II is a scatter plot of the optical and 
meteorological measurements of C^. The two lines drawn 
on either side of the equivalence line delineate the 

region where C^q equals within one standard deviation of 

experimental error. Sixty-five per cent of the data points 
are within this region indicating good correlation between 
measurements of C^q and C^. The meteorological values 
were systematically higher than the optical values. This • 
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would be expected in the case of saturation. However, the 
lack of a well-defined saturation level for optical values 
precludes a definite conclusion. 

This experiment indicated that the optical and 
meteorological systems that had been developed could 
measure atmospheric turbulence and, for the prevailing 
conditions of the experiment, the local thermal fluctuations 
formed a good predictor of beam behavior. 

B. SAN NICOLAS ISLAND EXPERIMENT 

The next experiment was conducted over the weekend of 
March 3 at San Nicolas Island. Laser beam scintillation 
measurements were taken over a one-kilometer path between 
San Nicolas Island and the R.V. Acania . This experiment 
tested the feasibility of making scintillation measurements 
between the shore and a ship. A helium-neon laser on the 
island was aimed at a shipboard detector package. Meteoro- 
logical measurements were made on-board the Acania . Inde- 
pendent measurements of C^q were made by the optical 
propagation group from the Naval Missile Center, Point Mugu. 

The data from this experiment are listed in Data Table 
II. Graph III shows the measurement of C^ as a function of 
time . 

The problem of keeping the laser beam lined up with the 
optical receiver proved very difficult. The laser operator 
aimed the beam at a retro-reflector mounted near the detec- 
tor and used the reflection for sighting corrections. 
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A gimballed mount was provided for the detector which could 
be aligned with the laser beam and thereby reduce the effects 
of ship motion on the detected beam intensity. This system 
proved crude and awkward to operate. By making measurements 
at night, the need for a background compensating photocell 
was eliminated. 

The major problem encountered was that the laser beam 
could not be pointed at the detector for periods of time 
longer than five to ten seconds due to ship motion. The 
detector was moved around in the beam by the roll and pitch 
of the ship causing a kind of inverse beam wander effect 
and theoretically increasing the measured values of C^q. 
Variations in laser beam intensity were due to both scintil- 
lation and ship’s motion. However, by keeping measurement 
times very short, the effect of ship motion was minimized 
and several log-intensity distributions were accumulated. 
Agreement to within about 40% of optical measurements along 
the beach line, optical measurements from shore to ship and 
meteorological measurements on shipboard lead to some 
confidence in the basic measurement techniques. 

This experiment showed, however, that a stabilized 
optical system was necessary to keep the laser beam and 
detector aligned.' The optical system, as constructed, could 
not reduce ship motion effects to an acceptable level. It 
was evident that the system needed to be re-engineered for 
shipboard use. 
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C. SECOND MONTEREY BAY EXPERIMENT 

The final propagation experiment was conducted on the 
second and third of May across Monterey Bay. As before, 
the laser radiation was transmitted from the Holiday Inn 
to Hopkins Marine Station over a horizontal path, five 
meters above the sea surface. The optical measurements 
were taken as previously described for the first Monterey Bay 
experiment. The meteorological measurements were taken on 
the beach at the Holiday Inn, where the temperature sensors 
were positioned approximately two meters below the optical 
propagation path. The values of were extrapolated for 
the difference in measurement levels. 

1 . Data Taken on May 2, 1973 

Administrative and technical problems caused many 
interruptions and delays in data accumulation during the first 
day of the experiment. Numerous beam alignment adjustments 
caused the system to be shut down often for realignment. 

The data collected, especially before 1800, were of ques- 
tionable accuracy. Graphs IV A and B show the optical and 
meteorological measurements of plotted as functions of 
time. The average values of the two graphs are approximate- 
ly equal, especially for data taken after 1800. The large 
variations in C^q are due primarily to beam alignment 
problems . 

2 . Data Taken on May 3, 1973 

The data taken on the second day of the experiment 
(Graphs IV C, D, E, F, G, H) showed better correlation 
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between the optical and meteorological estimates of C^, 
Between 0830 and 1030, the values of increased and 

varied rapidly in time. Corresponding measurements of 
Cj^q showed a lesser increase and did not reach the high 
values of C^. This was probably due to the saturation 
of Cj^q. Local cells of strong turbulence could have caused 
these effects but it was doubtful that isolated cells would 
have caused sequential increases in a period of two hours. 

Assuming that atmospheric turbulence over the 
propagation path was random, one expects that the number of 
localized cells of atmospheric turbulence that affect one 
point on the propagation path during a specific time period 
to be representative of the number of cells that passed 
.any point along the path during the same time period. 

There was no reason to expect that localized turbulence 
would cluster around the thermosonde system. One or two 
strong increases in C^ T during a long period of time could 
be blamed on isolated localized cells of high turbulence 
that would have little effect on laser beam scintillation. 
However, successive, strong increases in during a 

long period of time could only be due to an overall 
increase in atmospheric turbulence. These arguments tend 
to support the saturation explanation for the failure of 
Cj^jq to correlate with 

From 1200 to 1600, the values of C^q and 
showed high correlation, while atmospheric turbulence 
remained relatively steady. At approximately 1630, laser- 
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beam alignment was lost and data collection ceased until 
1845. During this time period, the values of showed a 

gradual decrease. When C^q measurements were re-commenced 
at 1845, values for C^q and C^ T showed a high degree of 
correlation. 

3 . Data Correlation 

Correlation of C^q with was checked by plotting 

Cj^q as a function of in Graph V (A, B, C) and Graph VI 

(A, B, C, D) . Graph V A shows the data measured on the 
second of May. Sixty-one per cent of all the data points 
lie within one standard deviation experimental error of the 
^NO'^NT ec l u i va l ence line. Experimental error bars have been 
plotted on only one point but apply for all points. No 
definite level of saturation was observed. Graphs V, B 
and C, show the data from the second of May sorted according 
to the mean wind speed at the time of measurement. No 
definite trends were noticed. For measurements made when 
the mean wind was less than seven miles per hour, 66% of 
all the data points lie within one standard deviation of 
the equivalence line. For measurements made when 

the mean wind was greater than seven miles per hour, 59% 
of all data points lie within one standard deviation of the 
Cno'^nt equivalence line. The estimated precision and the 
quality of the data and the lack of a wide range of wind 
speeds do not justify any significant conclusions concerning 
the effect of mean wind velocity on measurements. 
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Graph VI A is a scatter plot of the measurements 
made on the third of May. A large reduction in technical 
problems led to a higher level of confidence in the optical 
and meteorological data taken on this day. Sixty-three per 
cent of all data points lie within one standard deviation 
of the equivalence line. Graphs VI B, C, and D 

are scatter plots of measurements sorted according to the 
mean wind velocity at the time of measurement. For data 
taken during relatively calm periods (Graph VI B) , data 
point dispersion is very large. Only 231 of all the points 
lie within one standard deviation experimental error of the 
equivalence line. The data taken during periods when the 
wind speed was greater than three but less than six miles 
per hour (Graph VI, C) show much less dispersion. Seventy- 
nine per cent of the data points are located within one 
standard deviation experimental error of the equiva- 

lence line. This trend toward less dispersion’ with higher 
wind speed is continued in Graph VI D where data points 
measured during periods when the wind speed was greater 
than six miles per hour are plotted. In this case, 93% 
of all the data points are located within one standard 
deviation experimental error of the equivalence line. These 
data indicate that an increase in mean wind speed was 
accompanied by an increase in the homogeneity of the turbu- 
lence and thus, less dispersion of data points on the 
scatter plots. 
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VI. GENERAL RESULTS 



Empirical data from the Monterey Bay experiment indicate 
that the optical and micrometeorological measurements of 
atmospheric turbulence correlated well within experimental 
tolerances. Graph VII shows that 68% of all measurements of 
C^jq and were within one standard deviation experimental 
error of the equivalence line. The San Nicolas 

Island data were excluded from the statistical analysis 
because of the unresolved problem of ship motion. These 
results indicate that it is possible to predict laser beam 
scintillation at a wavelength of 0.63 micrometers in the 
marine boundary layer from measurements of mean pressure, 
mean temperature, and the covariance of the temperature 
fluctuations . 

Some evidence indicates that saturation of C XTri 

NO 

measurements may have occurred for values of C^j greater 

- 8 *1/3 

than approximately 5.0 x 10 meters . However, the 
lack of a definite level of saturation in the scatter plots 
suggests that definite conclusions on this must await 
additional measurements. 

The correlation between the optical and meteorological 
measurements of the refractive index structure constant 
increased as the mean wind speed increased. Measurements 
taken on the second of May were inconclusive due to a lack 
of range in wind speed values. It was expected that turbu- 
lence would become more homogeneous as the mean wind increased. 
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The theoretical model, from which the relations for 
were derived, assumed homogeneous turbulence in the hori- 
zontal plane. Hence, with an increase in wind speed, 
turbulence should become more homogeneous and the theoretical 
model should describe the physical processes more accurately 
in accordance with the experimental results. 

The San Nicolas Island experiment demonstrated the 
feasibility of making scintillation measurements between 
the R. V. Acania and a shore station. However, problems 
in laser alignment due to ship motion indicated that the 
optical receiver system needed to be redesigned. A 
stabilized control system is needed to keep the laser beam 
and detector system aligned. 

Finally, the high correlation between the optical and 
meteorological measurements of the refractive index structure 
constant indicates that Tatarski’s description of electro- 
magnetic propagation is essentially valid in the marine 
boundary layer at a wavelength of 0.63 micrometers. The 
validity of the theory appears to improve with increasing 
wind speeds. Additional experiments are needed to obtain 
conclusive evidence to support Tatarski’s formulation over 
a wide range of meteorological conditions. 
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VII. FUTURE RESEARCH SUGGESTIONS 



Refractive index structure constant measurements should 
be made over a wide range of turbulence to determine the 
possible existence and description of the saturation of 
optical values. These measurements could further substanti- 
ate or refute the existing results of this project. 

Project measurements were restricted to a single 
wavelength and propagation distance. Additional experiments 
should be made to determine the validity of the theoretical 
model's predicted dependency on the parameters of wavelength 
and propagation path length. Experiments could be conducted 
at various wavelengths using various types of lasers. 

Equipment should be designed to shield the laser from 
fluctuating meteorological conditions and to continuously 
monitor the laser's output intensity. This would reduce 
the internal intensity modulation and provide -a convenient 
method of determining beam stability. 

Additional thermosonde units should be added to the 
meteorological instrumentation to obtain a more representa- 
tive sampling of turbulence along the propagation path. 

This would reduce the possibility of localized isolated 
cells of inhomogeneous turbulence giving false indications 
of total path turbulence. Additional experiments should be 
conducted to determine if the shore line boundary has a 
significant effect on the meteorological measurements. 
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Simultaneous measurements of C xirp could be made both on the 
beach and on-board the Acania and compared for correlation. 

A digital system should be assembled for both optical 
and meteorological data. Digitizing the scintillation 
and temperature fluctuation signals would greatly reduce the 
experimental error and provide a direct interface with the 
Naval Postgraduate School's IBM computer. 

Finally, the development of a system to make ship-to- 
shore measurements of laser beam scintillation would be 
highly desirable. The system would require some form of 
stabilized optics to reduce ship motion effects. This 
would make the system totally portable and measurements 
could be made almost anywhere on the coast. 
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VIII. PROJECT CONCLUSIONS 



The experimental measurements obtained during this 
research project suggest the following conclusions: 

(1) Systems were developed that made possible simulta- 
neous optical and meteorological measurements of 
the refractive index structure constant in the 
marine boundary layer. 

(2) Laser beam scintillation at 0.63 micrometers was 
predicted from meteorological measurements of 
atmospheric temperature fluctuations, mean 
pressure and mean temperature. 

(3) Tatarski's formulation of electromagnetic 
propagation [1] was valid for propagation in 
the marine boundary layer at a wavelength of 
0.63 micrometers. 

(4) The validity of Tatarski's formulation increases 
with increasing mean wind speed and subsequent 
increasing homogeneity of atmospheric turbulence. 

(5) More research is needed over a wider range of 
meteorological conditions to substantiate the 
generality of the project's conclusions. These 
measurements would resolve the question of the 
possible existence and description of a satura- 
tion effect on optical scintillation measurements 
of the refractive index structure constant. 
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APPENDIX A 




Time Variation of the Refractive Index Structure Constant 



Measured with the Scintillometer System 




Time 10 January 1973 (PST) 

Time Variation of the Refractive Index Structure Constant •• 
Measured with the Thermosonde System 
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First Monterey Bay Experiment Data Correlation 



Graph II 




Refractive Index Structure Constant Measured with the 
Thermosonde System (X 10 ^m 

Data taken over Monterey Bay on January 10. Sixty-six per 
cent of all data points are within one standard deviation 



experimental error of the C 



NO NT 



equivalence line. 
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Graph III 
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0300 0400 0500 0600 0700 

Time 3 March IS 7 3 (PST) 

Measurements of the refractive index structure constant 
versus time. Data collected during the San Nicolas Island 
Experiment. 

° - measurements by NAVMISCEN optical propagation 

group across Laser Bay on San Nicolas Island. 

0 - Cj^q measurements made by NPS optical propagation 

group between, the R. V. Acania and San Nicolas Island. 

A - measurements by NPS meteorological group on- 

board the R. V. Acania. 
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Optical Measurements of versus Time of Day 




Time 2 May 1973 (PST) 

Meteorological Measurements of versus Time of Day. 
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Graph IV C 
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Time 3 May 1973 (PST) 



Optical Measurements of versus Time of Day. 




Time 3 May 1973 (PST) 

Meteorological Measurements of versus Time of Day. 
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Optical Measurements of versus Time of Day. 




Meteorological Measurements of versus Time of Day. 
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Graph IV G 
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Time 3 May 1973 (PST) 

Optical Measurements of versus Time of Day. 




Meteorological Measurements of versus Time of Day. 
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Graph V A 




Optical Measurements of C M versus Meteorological Measurements 
of C N . 



Data taken on the second of May over Monterey Bay. 



Sixty-one per cent of all data points are within one 
standard deviation experimental error of the 
equivalence line. 
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C NT (X 10" 8 m' 1/3 ) 

Optical Measurements of versus Meteorological 

Measurements of C XI . Data Taken when the Mean Wind 

N 

Speed was less than 7 MPH. 

Data taken on the second of May over Monterey Bay. 

Sixty-six per cent of all data points are within 
one standard deviation experimental error of the 
Cno" C nT ec ipi va l enc e line. 
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C NT (X 10 -8 m" 1/3 ) 

Optical Measurements o£ versus Meteorological 
Measurements of C^. Data Taken when the Mean Wind 
Speed was greater than 7 MPH. 

Data taken on the second of May over Monterey Bay. 

Fifty-nine per cent of all data points are within 
one standard deviation experimental error of the 
equivalence line . 



76 



10 



Graph VI A 




C NT (X 10' 8 m' 1/3 ) 

Optical Measurements o£ versus Meteorological 
Measurements of C^. 

Data taken on the third of May over Monterey Bay. 



Sixty-three per cent of all data points are within 
one standard deviation experimental error of the 
^NO'^NT ec i u i- va l ence line. 
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C NT (X 10' 8 m' 1/3 ) 



Optical Measurements of versus Meteorological 
Measurements of C^. Data taken when the mean wind 
speed was less than 3 MPH. 

Data taken over Monterey Bay on the third of May. 

Twenty- three per cent of all data points are within 
one standard deviation experimental error of the 
C N0" C NT ec l u i valence line. 



78 



10 



Graph VI C 




c NT (X 10'V 1/3 ) 



Optical Measurements of versus Meteorological 

Measurements of C... 

N 

Data taken when the mean wind speed was greater 
than 3 MPH but less than 6 MPH. 



Data taken over Monterey Bay on the third of May. 

Seventy-nine per cent of all data points are within 
one standard deviation experimental error of the 
equivalence line. 
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Graph VI D 
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C NT (X 10" 8 m' 1/3 ) 

Optical Measurements of versus Meteorological 
Measurements of C^. 




Data taken when the mean wind speed was greater than 
6 MPH. 



Data taken over Monterey Bay on the third of May. 

Ninety-three per cent of all data points are within 
one standard deviation experimental error of the 
Cnq-C^t equivalence line. 
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Graph VII 




C NT (X 10' 8 m" 1/3 ) 

Optical Measurements of versus Meteorological 
Measurements of C^. 

Includes all measurements of made over Monterey 
Bay. Excludes San Nicolas Island measurements. 

Sixty-eight per cent of all data points are within 
one standard deviation experimental error of the 
C N0~ C NT ec l u i va l ence line. 
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The Frequency Spectrum of the Scintillation Signal 
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Component Frequency (Hz) 



THE LASER SCINTILLOMETER SYSTEM 
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Figure 



LASER OUTPUT AND POINTING SYSTEM 



First lens 




This system was secured to a heavy metal table for 
vibration isolation. 



Figure 2. 
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OPTICAL DETECTOR PACKAGE 



0.63 Micrometer Bandpass Filter (30 X Bandwidth) 




To differential 
Amplifier 



Figure 3. 
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THE LASER SCINTILLOMETER 




Figure 4. 
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SCINTILLATION SPECTRUM ANALYZER 
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SCINTILLATION SIGNAL FREQUENCY COMPONENT INTENSITY VERSUS TIME 
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Figure 9. Laser Transmission Path over Monterey Bay. 




Figure 10. San Nicolas Island Experiment Transmission Paths. 
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APPENDIX C 



DATA TABLE I 

FIRST MONTEREY BAY EXPERIMENT 



C NO 


C NT 


— 


3.21 


2.87 


2.80 


2.50 


3.02 


— 


3. 56 


2.18 


2.51 


2.55 


2.80 


2.55 


3.22 


2.66 


3.60 


2.71 


3. 73 


2.98 


3.11 


3.30 


4.02 


3.57 


2.80 


4.36 


5.01 


4.15 


12.00 


— 


4.02 


3.25 


4.50 


3.41 


3.45 


— 


3.53 


3.25 


15.00 


— 


5.01 


— 


3.02 


3.00 


2.03 


2.87 


2.52 


2.93 


3.73 


2.77 


3.54 


— 


3.01 


2. 23 


3.70 


• 2.34 


3.00 


2.18 


3.55 
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DATA TABLE I 



FIRST MONTEREY BAY EXPERIMENT 



C NT 


C NT 


2.18 


3.11 


2.87 


3.52 


2.55 


2.91 


2.55 


2.43 


2.55 


2.72 


2.39 


2.84 


2.18 


2.45 


2.39 


2.61 


2.55 


4.10 


2.55 


3.04 


2. 55 


1.40 


2.61 


2.61 


2.61 


3.51 


— 


5.52 


— 


4.51 


2.66 


1.54 


2. 71 


3.40 


. 2.50 


2.70 


2.39 


2.51 


2. 23 


2.5 2 


2.29 


4.23 


2.13 


4.05 


2.18 


— 


. 2.07 


— 


2.39 


— 


2.07 


— 


2.02 


— 


2.07 


— 
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0645 



Time 

1400 

1405 

1410 

1415 

1420 

1425 

1430 

1435 

1440 

1445 

1450 

1455 

1500 

1505 

1510 

1515 

1520 

1525 

1530 

1535 

154.0 

1545 

1550 

1555 

1600 



DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



(2 May) 



'NO 



'NT 



— 


15.0 


— 


14.6 


6.01 


14.6 


— 


15.0 


5.60 


14.6 


— 


2.67 


5.52 


3.46 


5.60 


3. 46 


5. 79 


3.58 


5.60 


2.13 


5.01 


— 


5.34 


— 


5.25 


— 


6. 21 


::: 


- - - 


5.94 


— 


4.81 


— 


6.13 


— 


5.33 


— 


5.12 


5. 01 


2.67 


5. 69 


6.94 


5,90 


5.26 
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Time 

1605 

1610 

1615 

1620 

1625 

1630 

1635 

1640 

1645 

1650 

1655 

1700 

1705 

1710 

1715 

1720 

1725 

1730 

1735 

1740 

1745 

1750 

1755 

1800 



DATA TABLE III 



SECOND MONTEREY BAY 
(2 May) C NQ 



5.31 

7.01 

6.10 

6.72 

5.71 

6.25 

6.39 



5.70 



4.91 
4.07 
4.35 

4.92 



EXPERIMENT 

C NT 



5.33 

5.12 

4.61 

4.91 

4.29 
4.61 
4.61 
5.12 

4.30 
4.91 
4.81 
4.71 
5.12 
4.40 
5.02 
4.91 
4.81 
4.61 
4.74 
4.40 
4.71 
4.81 
4.30 
4.39 
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Time 

1805 

1810 

1815 

1820 

1825 

1830 

1835 

1840 

1845 

1850 

1855 

1900 

1905 

1910 

1915 

1920 

1925 

1930 

1935 

1940 

1945 

1950 

1955 

2000 



DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



(3 May) 



C 



NO 



5.35 

4.47 

4.77 

4.07 

4.49 
3.79 
3.51 
3.23 
3.65 
3.79 
4.57 

4.07 

4.07 

4.07 
3.93 
4.77 

4.50 

3.79 

3.65 

3.65 



C NT 



4.59 

4.59 

4.69 

4.59 

4.69 

4.69 

4.39 

4.59 

4.49 

4.39 

4.59 

4.69 

4.48 

4.28 

4.51 

4.48 

4.18 
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DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



Time (3 May) 


C NO 


C NT 


0630 


2.72 


— 


0635 


— 


— 


0640 


2.25 


— 


0645 


— 


— 


0650 


— 


— 


0655 


2.53 


— 


0700 


2.53 


— 


0705 


2.34 


— 


0710 


2.06 


— 


0715 


2.06 


— 


0720 


1.59 


— 


0725 


1.97 


— 


0730 


2.34 


2.31 


0735 


2.34 


3.08 


0740 


2.15 


1.95 


0745 


2.15 


_ 1.64 


0750 


2. 25 


3.60 


0755 


2.62 


3.60 


0800 


2. 62 


4.01 


0805 


2.28 


3.70 


0810 


2.15 


3.17 


0815 


2.62 


3.70 


0820 


2.81 


3.91 


0825 


2.72 


4.63 


0830 


2.72 


5.14 


0835 


3.47 


6.16 


0840 


3.56 


9.37 


0845 


3.56 


5.45 


0850 


3.56 


8.03 
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Time 

0855 

0900 

0905 

0910 

0915 

0920 

0925 

0930 

0935 

0940 

0945 

0950 

0955 

1000 

1005 

1010 

1015 

1020 

1025 

1030 

1035 

1040 

1045 

1050 

1055 

1100 



DATA TABLE III 



SECOND MONTEREY BAY 

(3 Ma Y) C N0 

3.84 

4.31 

4.12 

4.50 

4.31 

4.22 

3.84 

3.84 



2.35 

2.97 

2.15 

3. 09 
3. 09 
3.28 
3.28 
3.47 



4.22 

4.22 

3.71 

3.56 



EXPERIMENT 

r 

'-NT 

6.49 

5.81 

5.10 

6.12 
4.90 
7.97 

6.12 

6.12 

5.10 
4.39 

4.49 

5.10 
7.97 
5.61 

3.82 
4.25 
4.25 
8.85 
8.85 
4.70 
4.27 

3.50 
4.19 
3.84 
3.93 
3.42 
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Time 

1105 

1110 

1115 

1120 

1125 

1130 

1135 

1140 

1145 

1150 

1155 

1200 

1205 

1210 

1215 

1220 

1225 

1230 

1235 

1240 

1245 

1250 

1255 

1300 



DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



(3 May) 



C NO 


C NT 


3.47 


3.50 


3.47 


3.67 


3.47 


3.75 


5.15 


3.84 


5.90 


3.07 


— 


3.59 


4.12 


3.95 


4.03 


3.85 


4.12 


3.85 


4.87 


4.14 


4.78 


4.14 


4.31 


3.73 


4.22 


4.11 


3.94 


3.92 


4.40 


4.20 


4.03 


4.11 


4.69 


4. 39 


4.69 


4.20 


5.45 


3.83 


5.70 


4.48 


5.55 


4.67 


4.31 


4.67 


4.50 


4.47 


4.03 


4.29 
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Time 

1305 

1310 

1315 

1320 

1325 

1330 

1335 

1340 

1345 

1350 

1355 

1400 

1405 

1410 

1415 

1420 

1425 

1430 

1435 

1440 

1445 

1450 

1455 

1500 



DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



(3 May) 



C NO 

4.03 

4.03 

4.12 

4.31 

4.40 

4.59 

4.50 

4.78 

4.12 

3.65 

3.47 

3.47 

3.56 

3.56 

3.56 

3.56 

3.92 

3.75 

3.47 

3.56 



3.34 
3. 75 



C NT 

3.92 

3. 73 

4. 20 
4.29 
4.39 
4.50 
4.69 
4.41 
4.41 
4.69 
4.31 
4.50 
4.50 
4.50 
4.41 
4.78 
4.69 
4. 50 
4. 50 
4.60 
4.50 
4.22 
3.94 
4.22 
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Time 

1505 

1510 

1515 

1520 

1525 

1530 

1535 

1540 

1545 

1550 

1555 

1600 

1605 

1610 

1615 

1620 

1625 

1630 

1635 

1640 

1645 

1650 

1655 

1700 



DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



C3 May) C NQ 




4.12 

4.03 

4.03 

4.03 

5.25 



4.40 



5.15 

5.05 

4.69 

3.75 

4.40 

3.75 

4.12 



4.31 

3.94 

4.12 
3.75 
3.94 
4.41 
4.22 

4.13 
4 . 22 
4.22 
3.94 
4.12 
3.94 
4.31 
3.94 
3.75 
3.94 
4.24 
4.52 
2.82 
2.73 
3.18 
2.63 
2.82 
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Time 

1705 

1710 

1715 

1720 

1725 

1730 

1735 

1740 

1745 

1750 

1755 

1800 

1805 

1810 

1815 

1820 

1825 

1830 

1835 

1840 

1845 

1850 

1855 

1900 



DATA TABLE III 



SECOND MONTEREY BAY EXPERIMENT 



(3 May) 




1.97 

1.87 

1.78 

1.68 

1.87 



C NT 



3.01 
2.35 
2.07 
2.26 
2.35 

2.48 

2.48 

2.48 
1.88 
1.88 
2.06 
2.15 
2.24 
2.05 
1.86 
1.96 
1.86 
1.96 
1.68 
1.12 
1.68 
1.31 

2.01 

2.05 
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Time 

1905 

1910 

1915 

1920 

1925 

1930 

1935 

1940 

1945 

1950 

1955 

2000 

2005 

2010 

2015 

2020 

2025 

2030 

2035 

2040 

2045 

2050 

2055 

2100 



DATA TABLE III 



SECOND MONTEREY BAY 
(3 May) C NQ 



1.97 

2.06 

2.06 

2.25 

2.53 

2.72 

2.44 

2.53 

2.06 

2.25 

2.34 

2.34 

1.97 

1.78 

1.78 

1.68 

1.87 

2.06 

1.97 

1.87 

2.06 

1.96 

2.06 

1.69 



EXPERIMENT 

r 

^ NT 



2.34 

2.34 

1.96 

2.38 

2.01 

1.87 

1.68 

1.22 

1.40 

1.78 

1.80 

1.68 

1.73 

1.31 

1.78 

1.63 

1.78 

1.63 

1.78 

1.87 

2.10 

2.29 
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